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Abstract 



We discuss the implication of the introduction of an extra field to 
the dynamics of a scalar field conformally coupled to gravitation in a 
homogeneous isotropic spatially flat universe. We show that for some 
reasonable parameter values the dynamical effects are similar to those 
of our previous model with a single scalar field. Nevertheless for other 
parameter values new dynamical effects are obtained. 
PACS : 98.80.Cq, 98.80.Bp, 98.80.Hw 
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1 Introduction 



Recently, the authors investigated Q the dynamics of a spatially flat 
universe dominated by a self-interacting conformally coupled scalar 
field according to the action 

S=\J d^x^-g (- - + g^d^d^ - 2Vty) + W 2 ) , (1) 

where R denotes the scalar curvature, ip is the scalar fields, k = 8ttG 
(G being Newton's constant). A cosmological constant, if present, is 
incorporated in the scalar field potential 

vy>) = - —ip — 2 • ( 2 ) 

K 4 K 

Many new and interesting dynamical phenomena were discovered and 
its cosmological interpretation were discussed. Special attention was 
given to the conformal coupling case with £ = 1/6. This choice for the 
coupling constant is motivated by physical arguments from particle 
theories [Q] and from scale invariance at the classical level ||. The 
non-minimal coupling is also required by first loop corrections M. 
One interesting effect discussed in ^ is superinflation characterized 
by H > 0, and that can only be achieved by a non-minimal coupling. 

Here we consider the robustness of these dynamical phenomena 
with respect to the introduction into the model of a second massless 
scalar field. We also discuss some new type of asymptotic solutions 
arising from the introduction of the extra field. The structure of the 
paper is the following: in section 2 we present the model and some 
definitions. The fixed points and their stability are given in section 
3. Additional asymptotic solutions are discussed in section 4 and the 
paper is closed with some concluding remarks in section 5. 



2 The model 

We consider the conformally coupled theory described by the action: 
+ fdplhdvTk ~ 2W 2 ) + \b*& + Prfbl£\ ■ (3) 
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We use the full conserved scalar field stress-energy tensor 



V = d^ipidyijji — — (V^V;, — g^ u O) (tpf) + —G^tfjf 

+ d^ 2 d v ^ 2 - g (V^V^ - g^U) ($) + -G^l 

\g^{d^2d a ^ 2 -2V 2 {ij 2 ))-^\^l (4) 

(where G^ v is the Einstein tensor), thereby avoiding the use of any 
effective coupling constant in the Einstein equations. Moreover, this 
consideration of the energy-momentum tensor together with an 
adequate self-consistent treatment of Einstein's and Klein-Gordon equa- 
tions eliminates the artificial pathologies associated, in Einstein's equa- 
tions, with the "critical factor" 1 — (ipf + tp 2 ) K / The consensual 
attitude encountered in the literatures is indeed that the dynamics is 
ill-defined when this factor is negative or vanishes. As we will see, 
that this is not the case. 

Let us consider the following form for the potential V\ and V 2 : 



ViWi) = — fa - —^i - . (5) 



The parameter on is related to the mass of the particle by mi = 
yj 6a.i/n. As in the first papers, we are interested in the dynamics of a 
spatially flat Friedmann-Robertson- Walker universe with line element 
ds 2 = dr 2 — a 2 (r) (dx 2 + dy 2 + dz 2 ). Also we will restrict the potential 
in (H|) to the case a\ = a, a 2 = 0, Q\ = Q, £l 2 = 0, u\ = oj and uj 2 = 0. 
This corresponds to a massive field coupled to a massless field. The 
energy density and pressure associated to the scalar fields are: 

a = o x + cr 2 - -^lVl, ( 6 ) 

P = Pl+p 2 + ^1pl^ 2 , ( 7 ) 



where 



°i = Y + \ R2 ^ + \ H9 M) + vm, (8) 



Vl = & - \ [2HdMD + d 2 T M)] ~\{ 2 H + 2H 2 ) # - Vifa). (9) 
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The action (^) then implies the trace equation: 

R = -q(h + 2H 2 ) = — k (<j — 3p) , (10) 



the energy constraint: 

3H 2 -na = 0, (11) 
and the Klein-Gordon equations for the two scalar fields: 

$1 + 3HiPt - + ^ + PMl = 0, (12) 

^2 + 3Hijj 2 - \r^ 2 + ^ + 0fy 2 = 0. (13) 
6 aip2 

Using the Klein-Gordon equations (12) and ([[^) and the expres- 
sions for a and p we obtain 

a-3p = 4(V 1 + V 2 )-i> 1 f^- fa^. (14) 

dtpi dtp2 

This simple form holds for any interaction of the form /J^V^ provided 
7 + 5 = 4. Finally, the energy constraint is: 

3H 2 -ku = 3H 2 - ^Kipj - \^kB 2 ^\ - KiJVi^i - kVi 

~ \kH 2 iI% - kH^ 2 - kV 2 + = 0. (15) 

Now solving eqs. (0) and (|15|) for the derivat ives H and ipi we obtain: 

o 2 2 1 dVx 1 dV 2 

H = _2H 2 + -kV! + -kV 2 - ^i^f ~ r«^2TT. (16 
3 3 6 a^i 6 dip2 



ipi = - (-kH^! ± Vg) , (17) 



K 

where 

G = 6kH 2 - K 2 ip2 - 2k 2 H^2 ~ 2k 2 (Vi + V 2 ) - k 2 ^^ - K 2 H 2 ipl 

(18) 

The region G < in phase space is physically forbidden. 

Using eq. (17) allows to rewrite system (p|-|l~0|) in the following 
equivalent form: 

H = -2H 2 + ^-^, (19) 
b 



^ = -H^ ± (20) 

$2 = -3H^ 2 + \r^ 2 - ^ - W^2- (21) 
o aip2 

The system of ODE (ffap|gD is a five-dimensional system, but due 
to energy constraint (|ll[) the real phase space is constrained into a 
four-dimensional manifold. The situation is similar to the one field 
case, where the original system is three dimensional, but the energy 
constraint enforces the motion to take place on a two dimensional 
manifold. In the one field case, the plane (H, ip) is divided in sectors 
by some straight lines, corresponding to regions of distinct state equa- 
tions for the fluid ip. In the present case, the projection of the real 
phase-space onto the (H, ipi) plane conserves some of that sectors. For 
example, from 

a - 3p = 6a*l%, (22) 

valid for any quartic interaction potential between tpi and ip2-> one can 
see that the region ipi = still corresponds to a state equation of pure 
radiation. 



3 Fixed points 



In our case De Sitter solutions correspond to fixed points except for 
the fixed point at the origin. The fixed points are obtained from the 
conditions ip\ = ip2 = 0. The stability of the fixed points is determined 
from the linearized equations in a neighborhood of each point, except 
for special cases discussed below. For all the fixed points we have, 
in the formula below, si,S2,S3 = ±1. The existence conditions for 
the fixed points is such that the arguments of the square roots are 
positive: 



Type A 




(23) 



/6 Mfl+jMajJ 

^ 2 = S3 VkV WTW) ■ (25) 
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• Type B 



, (26) 

K 

ih = i>2 = 0. (27) 



F = Sl J^ A /__^. (28) 



a 2 - 







— a 


/ a 


— a; 




— a 


0. 





• Type C 



V>2 = 0. (30) 
Fixed points of type B are on the G = surface. 

In order to discuss the stability of these fixed points let us consider 
the generic case given by a system of the form 

Xi = Fi(xi, . . . ,x n ) = F(x), i = l,...,n. (31) 

A fixed point x satisfies Fi(x) = 0. It stability is usually established by 
linearizing ( |3l[) in a neighborhood of x. In this way we write x = x+5x 
and obtain (|3l|): 

'' l):r J5x + 0(5x 2 ), (32) 



dt 

where J is defined by 

ni- 

(33) 



J = ®* 

lJ dxj 



and Sx = (Sx±, . . . ,5x n ). If Det(J) ^ and if the eigenvalues Aj of 
J are such that Re(Aj) ^ then the Hartman-Grobmann theorem 
ensures that the linearized system obtained by retaining only linear 
terms in ( |32"| ) is topologically equivalent to the original system ( |3l] ) in 
a neighborhood of the fixed point. The real parts of the eigenvalues 
then determine the local stability of the fixed points. The stability 
of the fixed points on the G = surface cannot be studied using this 
method as Det(J) = 0. To overcome this problem we will consider 
the system formed by eqs. flT2"D, fll3| ) and (|l6|). The only points to be 
considered separately are of type B and their stability is studied by 
other methods. 

Equations (|l^) and ( |l3|) can be rewritten as a first order equations 
by introducing the new variables <p\ = ip\ and <p2 = ip2- In this way 
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we identify the variables X{ in ( pip with (H, ipi,4>i,ip2,(l>2) and obtain 
from (j33|) the following form for the J matrix computed at the fixed 
point of coordinates (H, tjj 1 , fi-y, ifi^, (j> 2 ): 



(34) 



/ 


-AH 2o^i 








o \ 







1 










Ci 


-3# 


















1 


V 


-2{a + P)^ 2 





Ci 


-3H J 



with 



and 



Ci 



-Pijj 2 2 + 3(0 - a)ipl + -{oj- a), 

K 



Co 



-{a + p)$\ + 6-. 

K 



The eigenvalues of J are: 
A = 



-H±-\ 9H + 9A±-VB, 
2 2V k ' 



A 



(35) 
(36) 

(37) 
(38) 



where the ± signs are to be considered independently of each other, 
A = |(3fi - j3 - Aa)^\ + 4a)S + ^(2w - a), (39) 

and 



= k 2 (/3 2 ^2 + (20a/3 + 14/3 2 - 6/30)^2 
+(4a 2 + /? 2 - 4a/3 - 120a + 6/30 + 90 2 )^) 



(40) 



A fixed point is stable, unstable or a saddle point if the real parts of 
the five eigenvalues are all negative, all positive or have different signs, 
respectively. Points of Type A and C are saddle points. The stability 
of fixed points of type B is discussed in the next section 



4 Asymptotic dynamics 

Two kinds of asymptotic dynamics are more relevant for our model: 
attractive fixed points and diverging solutions. For the one field model 
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the fixed point at H = ip = acts as a an attractor and plays an 
important role in our approach ||]. The fixed point at 

ipi = fpi = 1P2 = 1P2 = and H = a/— 3u>/k, (41) 

has a similar role. For this purpose we will consider here the special 
case of a non-negative cosmological constant u) < 0. This implies that 
H > for H = as a — 3p is positive for u negative and therefore the 
region H > is invariant under the dynamics. Let us now consider 
the function 

l = U\ + \^-i>\ + + -M - n)^ - \mi (42) 

z zp-a k 4 z 

Its total time derivative modulo eqs. (|l2j), (|l3|) and ([!(]) is given by 



L = -3F + tt^V>2 + — ^ + x ^2- (43) 

\ p — a J k (p — ajK 

For parameter values such that u = 0, (3 < and < a the func- 
tion L is a Lyapunov function for the fixed point ( f4l| ) and satisfy the 
properties 

L > 0, (44) 

L < 0. (45) 

The equality in (^4|) and ( |45| ) is valid only on the fixed point. Hence 
the fixed point is a global attractor for the whole region H > (see 
fig. 0). Numerical simulations indicate that the same is also true for 
uj < as shown in fig. |2|. 

For other parameter values numerical solutions point for an asymp- 
totic diverging solution (see fig. ||). For simplicity, let us consider the 
case uj = and the ansatz: 



H = 5i/>i, (46) 



with 6 a constant. Hence 



H = 54>i. (47) 
From eqs. (|l6|), ([l?]) and ( f47| ) we obtain 

<tyi = ^?(-12£ 2 + m? K ), (48) 
o 

which is a closed equation for ipi, with solution 

^ {t) = 65 + (m^-HnDM^t ' (49) 
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la lb 




Figure 1: Example of solutions projected in the (H,i[>i) plane for u = 0. 
The parameter values used are a = 1, Q = 1/2, m = 1, (3 = —1.2345. The 
solutions are attracted to the fixed point at the origin. 



2a 2a 




Figure 2: Example of solutions projected in the (H,ipi) plane for uo = —0.5. 
The other parameter values are the same as in fig. |l|. The solutions are 
attracted to the fixed point at iftx — and H = a/ —3uj/k. 
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3a 3b 




Figure 3: Solution with u = —0.3. The remaining parameters are as in fig. |l|. 
Some solutions are diverging while ore are still attracted to the fixed point 
at the origin. The quantities H, tpi and i[)2 diverge in a finite time. 



where ipi(0) is the initial condition for ipx at t = 0. Using eq. ( |12| ) and 
solution (H) yields easily ip2{t)'- 



. . , l~2y/365 A - 125 2 m(K - 18Q5 2 + mf K 2 z/>i(0) 
^ 2(t) - ± V"^ (1252-«m?)^(0)t + 65 " (5 ° } 

With the expressions for V'i(^) and it is straightforward (but 
a little bit cumbersome) to obtain the proportionality constant 5 from 
the remaining equation (0): 



§ = ±iy6/em? -9/3 + 3^-12/3^ +9/? 2 . (51) 

This value for 5 agrees with the results of numerical integrations. It is 
important to note at this point that not necessarily every solution has 
the asymptotic behavior of eq. ((!(]). From numerical investigations 
it seems to exist a threshold in the initial values of ip2 and ip2 such 
that above it the solutions converge asymptotically to the behavior 
obtained above. Below this threshold solutions near the origin in the 
(H,ipi) plane are attracted towards it, exhibiting a similar behavior 
to the case of a single massive scalar field. 



5 Concluding remarks 

The introduction of an extra massless field conserve some essential 
features of the dynamics of the one- field model of ref. Q, for some 
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parameter values, particularly the spiraling solution near the fixed 
point at the origin. Nevertheless some new type of solutions exist that 
have no analog in the one- field case. For two self- interacting coupled 
fields the behavior is even more richer and a more exhaustive study 
is under preparation. Nevertheless, for some parameter ranges, new 
dynamical effects are obtained which modify the model in a substantial 
way. A discussion of which parameter values are of physical relevance 
is important for the future development of the present approach. 
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